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We develop a microscopic toy model for Cavity dynamical Casimir effect (DCE), namely, the 
photon generation from vacuum due to a nonstationary dielectric slab in a fixed single mode cavity. 
We represent the slab by A 3> 1 noninteracting two-level atoms coupled to the field via the standard 
dipole interaction. We show that the DCE is contained implicitly in the light-matter interaction 
Hamiltonian when its parameters are externally prescribed functions of time. We also predict sev¬ 
eral new phenomena, such as saturation of the photon growth due to effective Kerr nonlinearity, 
generation of pairs of atomic excitations instead of photons (“Inverse DCE”) and coherent annihi¬ 
lation of pair of system excitations due to the atomic modulation (“Anti-DCE”). These results are 
extended to the circuit QED architecture, where similar effects can be implemented with a single 
qubit providing an alternative way to generate cavity and atom-held entangled states. 
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I. INTRODUCTION 


The term “dynamical Casimir effect” (DCE) is used nowadays for a rather wide group of phenomena characterized 
by creation of quanta from the initial vacuum state of some field due to time-dependent variations of the geometry or 
material properties of a macroscopic or mesoscopic system (see [Il-@ for recent reviews). In particular. Cavity DCE @ 
denotes the process of photon generation from the electromagnetic vacuum (and other initial states) in cavities due to 
the motion of some wall or the time-modulation of the material properties (e.g., dielectric permittivity or conductivity) 
of the wall or a medium inside the cavity (3)0- An analog of Cavity DCE was recently iimlemented experimentally 
in the solid state architecture known as circuit Quantum Electrodynamics (circuit QED |lMlll|). where a Josephson 
metamaterial was embedded in a low-Q microwave cavity, permitting the modulation of the cavity effective length 
via external magnetic field [l3 |. 

Although Cavity DCE has been studied theoretically for more than four decades, some aspects of this phenomenon 
are still not completely clear. A particular issue we approach here is the asym pto tic behavior of photon generation: 
while some models predict the saturation of the intra-cavity photon number [131 Il4| , other predict exponential photon 
growth even in the presence of moderate dissipation [1, [isMT^ . This controversy can be resolved by constructing a 
full microscopic model for the interaction between the quantized electromagnetic field and moving or time-modulated 
objects constituted of individual atoms. Some steps along this line were taken in [mu, yet the majority of studies 
employs time-varyin g bo undary conditions for the cavity field to bypass the complicated light-matter interaction at 
the interface [l], jj, □, I22 h^ . 

In this paper we utilize the general mathematical description of nonstationary circuit QED systems formulated 
recently in [28| to develop a microscopic toy model for Cavity DCE. Our study is motivated by the following intuition: 
since the boundary conditions are just a mathematical artifact to manage the interaction between photons and 
a large number of atoms, DCE should ultimately originate from the basic form of light-matter interaction with 
nonstationary parameters. So we consider the special case of Cavity DCE implemented with a dielectric slab having 
externally prescribed motion and material properties. The slab is portrayed as an ensemble of N two-level atoms 
with unspecified transition frequencies and coupling strengths that interact with the field via the standard dipole 
Hamiltonian with time-dependent parameters [29|. We use the time-independent boundary conditions to quantize 
the cavity field in a standard manner, while the interaction between the arbitrarily modulated atoms and photons is 
treated microscopically. 

After cumbersome calculations we arrive at simple mathematical expressions that generalize the common DCE 
description in single-mode cavities Hill. In particular, we express the photon generation rate in terms of the 
microscopic parameters, show that the photon growth and amount of squeezing are limited due to effective Kerr 
nonlinearity and point out that Cavity DCE occurs even for a single atom. Moreover, we discuss how external 
classic al p umping can significantly enhance the photon generation from vacuum for suitable choices of the pump 
phase 3]J. Since our model is quite genera,!, we also apply it to situations where all the system parameters are known, 
such as a cloud of cold polar molecules or superconducting qubits [Tl|,[ll,[3^. New effects arising from periodic 

external modulations are analyzed: generation of pairs of atomic excitations from vacuum (“Inverse DCE”), coherent 
annihilation of a pair of system excitations (“Anti-DCE”) and generation of entangled light-matter states. 
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This paper is organized as follows. In section |n] we formulate our problem and in section IIIII we develop the toy 
model for cavity DCE, presenting the analytical and numerical results. In section HVl we extend our analysis to cold 
atomic clouds, where all the atomic parameters are known and, in principle, can be modulated externally. In section 
El we repeat this analysis for the case of a single two-level atom, discussing the Anti-DCE behavior and studying 
some applications in the area of circuit QED. Section IVTl contains the conclusions. This paper contains two extensive 
appendices: in|A]we give the thorough analytical description for the case ^ I in the Heisenberg picture, while in 
m we do the same for iV = I in the Schrodinger picture. 


II. MATHEMATICAL FORMULATION OF THE PROBLEM 

We quantize the cavity field using the standard methods with time-independent boundary conditions [nisi]. The 
annihilation and creation operators a and do not depend explicitly on time, so the vacuum state defined as a|0) = 0 
0,i is the same for all times, unlike the case of a cavity with moving walls for which the field state depends on the 
instantaneous frequency (l^ . We consider a small dielectric slab located at an arbitrary position within the cavity, as 
depicted in figure I. The dielectric slab is subject to pre-determined motion with small amplitude, and its material 
properties (e.g., dielectric permittivity) can be adjusted externally by some bias (represented by the laser beam in the 
figure). From the microscopic point of view, this problem corresponds to a predetermined motion of an atomic cloud 
whose internal properties are modulated externally. For consistency, the generation of photons from vacuum in this 
particular example of DCE should be contained intrinsically within any formulation of the light-matter interaction. 



Zoom 


FIG. 1: Artistic view of DCE due to a nonstationary dielectric slab in a fixed single-mode cavity. The dielectric slab (pictured 
as a set of N noninteracting Hydrogen atoms) oscillates according to an external law of motion, while its dielectric properties 
can be modulated externally via electric or magnetic fields. The harmonic wave represents the time-independent cavity mode 
fnnction; red beam represents the modulation of the material properties of the dielectric. The zoom shows an individual atom 
containing one proton and one electron, whose center-of-mass coordinate R changes due to the prescribed motion. 


We consider the simplest microscopic model for the dielectric slab - a set of N non-interacting Hydrogen atoms, 
as shown in the zoom of figure I. First we recapitulate the interaction of a single atom with the field. Each atom 
consists of a proton (electron), described by the position operator fp (fg), with mass nip (mg) and charge e (—e). 
Introducing the center-of-mass (CM) position operator R = {nie^e + 'mprp)/M, where M = rUg -I- nip is the total 
atomic mass, we define the momentum operator P = Pg -|- Pp associated with the CM motion, where Pp (pe) is the 
canonical momentum operator of the proton (electron). Furthermore, one introduces the relative coordinate between 
the proton and electron f = fg — fp and the momentum p = (wpPg -I- mePp)/M associated with the relative motion 
of the reduced mass fi = nienip/M. As a result, we can decompose the dynamics into the motion of CM and the 
relative motion, with the total kinetic energy given by Pg/(2me) -I- Pp/(2TOp) = P^/(2M) -|- p^/(2^). 

We treat the light-matter interaction in the first-order dipole approximation, assuming that the dimensions of the 
atom are much smaller than the wavelength of the cavity mode. The minimal coupling Hamiltonian in the Coulomb 
gauge is minutely deduced in [2^. Considering that the CM motion is prescribed externally, with P and R given by 
known functions of time, it reads 


hW 


2M 


+ Hf + Ha - er ■ 


E -f R X B 




(f.B)^' 
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e(wp - We) 

2M 


■ B • (f X p) + (f • VR)f • E 


( 1 ) 


Here Hf = hujfi is the cavity free Hamiltonian, where uj is the frequency and n = a)d is the photon number operator. 
Ha = P^/(2m) + ^(f) is the Hamiltonian of the atomic internal dynamics, where V (f) = —e^/(47r£o|f |) is the Coulomb 
interaction energy and Eq is the permittivity of vacuum. The electric and magnetic intracavity fields are 


E (R) = z 



(R) (a — fii) 


( 2 ) 


B(R) = y^^(Vxu(R))(a + at) , (3) 

where V is the mode volume and u (R) is the dimensionless mode function determined from the time-independent 
boundary conditions on the walls. In the stationary case, when R = 0, one usually neglects the contributions con¬ 
taining the magnetic field and the gradient of the electric field in Hamiltonian ([IJ, recovering the standard dipole 
interaction term —ef • E. However, in the nonstationary regime R 0 all the terms must be taken into account [^ . 

For our toy model we take into consideration only the two atomic levels near-resonant with the cavity frequency, 
restricting the atomic dynamics to the “ground” and “excited” states \g) and |e), respectively. Hence the atomic 
Hamiltonian reads Ha = where is the transition frequency and = |e)(e| — \g){g\ is the Pauli operator. 

For the two-level approximation to hold we must have |w — f2| <C w. The position operator can be written as 
f = ro(T+ d-rpo--, where rg = (e|f I5) is the off-diagonal matrix element and the Pauli ladder operators are = \e){g\ 
and (7_ = \g){e\. In this case = |rop and the square of the magnetic held operator oc (a -I- appears naturally 
in Hamiltonian O- 

Hence the simplest model for a nonstationary dielectric slab in a stationary cavity is described by the general 
Hamiltonian of the form (we set ?i = 1) 


N r 

H = Lon -f 

1=1 


n 


.(0 


■5(c 


xt)(d« 


+ d«) 


-I- ~ + id{d^ — d) , 


( 4 ) 


where the index I labels the identical noninteracting atoms. The renormalized cavity frequency uj is constant, while 
the atomic transition frequency H, the atom-cavity coupling strength g and the “squeezing coefficient” y are regarded 
as externally prescribed functions of time. This occurs both due to the motion of the slab and the external in situ 
modulation of the atomic properties, though here we do not pursue the exact dependence. The last term on the 
right-hand side (RHS) of equation (|1]) accounts for the classical one-photon pumping of the cavity field @, included 
for generality and to study how DCE can be enhanced by an additional coherent drive. 

To understand the emergence of DCE from the microscopic viewpoint we do not need to know the exact relation 
between the parameters of Hamiltonians (ED and o, since for a weak external perturbation of the system we can 
write 


X = Xo+exY^ sin -k , ^ = {w, D, g, y, d} , (5) 

where Xq is the bare value and ex > 0 is the modulation depth of X. The sum runs over all the present modulation 
frequencies we can write it as +Y''j^ where Y'j denotes the sum over “fast” modulation frequencies, 

> Wo, and Y'j ~ “slow” modulation frequencies, <C wq. Parameters > 0 and the relative 

weights and phase constants corresponding to the modulation of X at frequency g^^^. For the classical pump we set 
do = 0, and we included the modulation of w in equation (O for the sake of generality. For the future use we define 
the complex modulation depth that includes the weight and the phase of X-modulation at frequency g^d) 

= £xw%^ exp[#^^] , X = {uj, D, g, y, d} . (6) 

Throughout the paper the notation and stands for the complex modulation depths corresponding to fast 
and slow modulation frequencies, respectively. 
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III. TOY MODEL FOR DCE WITH A DIELECTRIC SLAB 


As a toy model for DCE we consider a fixed cavity of known frequency ujq that contains N identical two-level atoms. 
The atomic transition frequency D and the coupling strength g are unknown, but in order to represent the dielectric 
slab the difference Wq — D must be large compared to the coupling strength, jwo — D| l^j. Due to the external 
perturbation the parameters D, g and y vary according to equation ([5]), and we consider the general case when the 
three parameters can change simultaneously. For the macroscopic slab we consider N ^ 1 and define the collective 
operators via the Holstein-Primakoff transformation [s^ 

=b\N = {N (7) 

1=1 1=1 1=1 


where the ladder operators h and w satisfy the bosonic commutation relation [&, 6 ^] = 1. To the first order in wh/N 
the Hamiltonian for our toy model reads 


H = ijJqu + -|- g{a + a^){b + 6 ^) + — a^) 


^ia + a^){bH + Pb^) , 


( 8 ) 


where we dehned the collective coupling constant g = i/Ng, so that go = \fNgo and Eg = yfNsg (we consider 50 > 0 
without loss of generality). In this paper the tilde over a c-number corresponds to the collective A-atoms parameter. 
The Hamiltonian (|S]) holds provided the inequality {wb) <C A is satisfied. 

In the dispersive regime, |A_|/2 ^ go, where A_ = wq — Dq is the bare atom-field detuning, the approximate 
solution in the Heisenberg picture is deduced in I A 21 


a 

h 


„-i(wo-<5+-<5x)* 




A_ 


B - sin 

A_ 



( 9 ) 

( 10 ) 


A_|_ — ojq -|- Dq 


< 5 + = 


9o 


ojQ ± Do 


— 


4Xo 

Wo + Do 


( 11 ) 


A and B are independent bosonic ladder operators that obey the Heisenberg equation of motion idO/dt = [O, Ag//] 
(O = A, B) with the effective Hamiltonian 


^eff = He + Hng ■ 


( 12 ) 


Here Hq contains the Gaussian part (quadratic terms in the operators A and B) and Hj^q contains the non-Gaussian 
part (quartic terms). 

For the modulation frequency 


= 2(wo -I- - 5+ - (5^ - C) 


(13) 


where we introduced the small adjustable “resonance shift” C in order to perform the fine tuning of the modulation 
frequency, we find 




Ho = qe 


^2 AB -f f — ) 

A_ V A 


h.c. 


(14) 


q = 


iS-Qo 

(wo + Do) 


-iD) ~{D) 


2Do go 


-(D) 


(15) 


Hng = — 




Mit 2 i 2 _ ^^ 2^2 ^ 2iliBlB 

A?. A_ 


(16) 
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^9o Q2^-2iS^t _ 4g0 J\i2 

A'i A_ 


+ i 


.5. e'^) 


2N go 


g g- 2 *ct 


A^B^B 


+ ^i1'i3 + - -p-A^B^ 

2Al_ A_ A_ A_ 


A 2 


~2 ~3 

-t _j__^0_^t^3g—_ _^2_y^t^3g —3i5^t 

A^ A^ 


/l.C. 


These results were obtained under a series of approximations. First, the detuning and the modulation depth of the 
atomic transition frequency must be small, £n, |^-| ^ wq, while the modulation depth of the atom-field coupling 
strength is ig go- Second, there are some restraints on the number of excitations in the atoms-field system for 
which our approach is accurate: 


N 


\J{(A a) -I- (&tS) 


go gg Ixol £x £o£x 1 , 

Wo ’ A_ ’ Wo ’ Wo ’ A_wo ’ A?_ J 


(17) 


^\l (a^a) {Wh) 


(jAh) \J go (a’^a) 
A ’ N ’ |A_| N 


< 1 . 


(18) 


Third, in equation (IT^ there are small “Systematic error frequency shifts” (SEFS) Ag that were neglected in order 
to keep the formulae concise. They are of the order 



Hence in the actual implementation of DCE one has to find experimentally the exact modulation frequency by scanning 
within the range A? 7 , so in part we introduced the adjustable resonance shift C to achieve this fine tuning. 

One can simplify the Hamiltonian (fT^ a little further. Neglecting the non-Gaussian terms we have 



go g 
A_ dt 


( 20 ) 


so for Ig| 1(5^1 we can write 


B{t)-B{0) + e^*^^^ 


Ait) - i(0) 


( 21 ) 


Assuming that the cavity and the atoms were initially in the ground states and substituting equation m into m, 
to the lowest order in go/A_ we obtain the Hamiltonian 


H, 


eff ^ qe 


-h /l.C.) - Na f(iU)2 -h iU 


( 22 ) 


where 



(23) 


is the effective Kerr nonlinearity strength due to a single two-level atom. Defining the phase (j)q via the relation 
q = i\q\ and introducing the new annihilation operator 

dr = (24) 

(that also satisfies [dr, d[(] = 1), the evolution of dr is governed by the time-independent Nonlinear DCE Hamiltonian 

HdCE = UJrfLr + arfll + igr(d|^^ - d^) . (25) 
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Here fir = ajar, Wr = (C — Nee), ^ = —Na and qr = Ig]. The term igr(dj^ — aj), which describes the simplest 
case of DCE in oscillating cavities Q, appears naturally in our derivation. The Hamiltonian (1^51) is well known from 
Nonlinear Quantum Optics for describing (in the interaction picture) a cavity that contains a Kerr medium and is 
parametrically driven , so we call ujr “effective detuning”. 

Hence we were able to deduce microscopically the DCE from the most basic form of light-matter interaction, 
equation ([8]). It turns out that DCE is described by the cumbersome non-Caussian Hamiltonian given by equations 
(fTl) - (USD, and the standard expression for cavity DCE is recovered only to the lowest order in go/A^. Recalling that 
the auxiliary annihilation operators A and B are related to the physical annihilation operators a and b via relations (|^ 
we can formulate the first new prediction of our toy model: the photon creation from vacuum is accompanied 
by the excitation of the internal degrees of freedom of the atoms in the slab, which becomes entangled with the cavity 
field. As stated previously, we assume that the CM motion of the atoms is prescribed externally, so our model does 
not contemplate the important back-action effects of DCE on the motion of the slab ilEaiEIl. 

The simplest realistic description of Cavity DCE must include (at least) the Kerr nonlinearity, as shown by equation 
(1^51) . Although separately the DCE and Kerr Hamiltonians can be integrated in a straightforward manner [^, the 
general analytical solution for the nonlinear DCE Hamiltonian is not known. To get qualitative insights about 
the asymptotic dynamics of Hamiltonian (1251) we rewrite it in the form of interaction picture parametric amplifier 
Hdce = Drfi-r + iqr{al^ — Q-J), where the overall detuning operator is Dr = uJr + arfir- Treating the detuning as a 
c-number {Dr), the solution in the Heisenberg picture reads 

ar {t) = T*ar (0) Gal (0) (26) 

D = cosh (Bt) + i sinh (Bt) , C/= 2^ sinh (Ht) , B = \J 4gJ — (Dr)^ ■ (27) 

So we arrive at the second new prediction of our model: asymptotic exponential photon growth is impossible for 
nonzero Ur, as for any fixed value of LOr the parameter B becomes imaginary for {fir) —> oo. This results solves the 
controversy about the long-time behavior of Cavity DCE, supporting the finding [IJ, [IJ that the photon generation 
is limited even in the absence of dissipation. 

To elucidate the system behavior for finite {hr) we write the wavefunction corresponding to the Hamiltonian Hdce 
as 


\ip) = eyi'p[—it{ujrrn + arrn^)\crn\rri) , (28) 

m—0 

where \m) denotes the Eock state. The probability amplitudes obey the differential equation 

Cm = g,[v'm(m- - sj{m + l)(m + 2)e-2*‘[“^+"“^(’”+i)lc„,+2]. (29) 

One can easily solve the pair of equations connecting just the amplitudes Ck it) and ck +2 (t) [1^. Eor ck +2 (0) = 0 
we get 


CK+2{t) = qrViK+ l)iK^ 

Hk 

Rk = v'K + 2a,(K + 1)]2 + g2(if + i)(if + 2) . (31) 

So the probability amplitude ck +2 is approximately decoupled form ck when qre/(K + l)iK + 2) <C Rk- For 
2\ar\ {K + 1) Iwrl the decoupling condition becomes qr <C 2|q!i.|. Therefore, in order to generate many photons 
from vacuum one must satisfy the condition qr ^ 2|ar|- In the opposite case, qr ^ {cerl, we expect generation of a 
small number of photons [l^. For example, for AT = 0 and qr/\ar\ <C 1, the effective detuning must be adjusted to 
LOr = —‘2ar to Optimize the coupling between the probability amplitudes cq and C 2 . As qr/\ar\ increases one can set 
LOr = —2 {K -I- 1) ar to Optimize the coupling between the amplitudes {ck, ck ^2 \ {K = 2,4,...), while the off-resonant 
coupling between {cq, C 2 , • • •, Cic} still allows a substantial population of ck [53|. So the question of utmost practical 
interest is: what value of tOr, or equivalently, what value of the adjustable resonance shift C optimizes the photon 
generation from vacuum in the presence of Kerr nonlinearity? The answer will given in the next subsection with the 
help of numerical simulations. 
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A. Numerical results 


We studied numerically how the Kerr nonlinearity affects the photon generation from vacuum. For the sake of 
completeness we included the cavity damping by means of the standard master equation at zero temperature 


dp 

dt 


-i[HDCE, P] + ^ {2drpdl - dldrp - pdldr) , 


(32) 


where p is the density operator, k is the cavity damping rate and Hdce is given by equation (El]). Strictly speaking, 
the microscopic derivation of this master equation does not contemplate the nonstationary case studied here, when 
the system parameters vary rapidly with time and the counter-rotating terms play a fundamental role (53| . Hence the 
solution of the master equation can only be used to grasp qualitatively the overall effect of dissipation. The stationary 
state of equation (15^ can be calculated analytically using the method of potential solutions for the corresponding 
Fokker-Planck equation |4^. However, as shown in figure [2l the asymptotic solution is of little help for our 
problem because the cavity field state during the time period of interest (initial times) may be very different from the 
asymptotic one. 



FIG. 2: Time behavior of the average photon number (fir), Mandel Q-factor and the variance of the squeezed quadrature 
(Apr)^ obtained via numerical integration of equation (I32II for qr = 3|Q:r|. For curves 1-5 the initial state is the vacuum state. 
For K = 0 the curves are: cJr = 0 (1), u)r = —Sor (2), u)r = —lOor (3) and tUr = —12Q:r (4). Line 5: uir = —lOor and k — qr', 
the dot-dashed line indicates the asymptotic value in this case. The dashed line corresponds to the initial thermal state with 
the average photon number n = 0.1 and parameters ujr = — lOor, k = 0. 

For small ratio qr/\ar\ “C 1 only two photons are generated from vacuum for ujr = —2ar, as predicted by equation 
(1501) . For larger ratios qrl\ar \ the behavior becomes much more complicated and the dynamics is shown in figures [5] 
and [5] for different values of the effective detuning ujr (which can be adjusted experimentally by tuning the resonance 
shift C). We plot the time behavior of the average photon number {hr), the Mandel Q-factor and the variance of the 
squeezed field quadrature (Ap^)^ = (p^) — (Pr)^, where 

{hr{hr - 1)) - {hr)"^ . 0^-4 

Q = -- ) Pt = —— . (33) 

{Ur) 2,1 

In figure [5]we set qr = 3|ar| and in figureOgr = 650|ar|. As expected, by increasing the ratio qrl\ar\ more photons 
are created from vacuum, and {hr) can be optimized by choosing an appropriate value of uJr- The average photon 
number is limited from above and exhibits a sort of irregular collapse-revival behavior due to the Kerr nonlinearity, 
as opposed to the exponential photon growth for the pure DCE case [1, ITsl - fT^ . The quantities Q and (Apr)^ also 
undergo oscillations, but they do not collapse to their initial values, meaning that the field state never returns to the 
vacuum state. The collapse-revival behavior of {hr) was discovered more than two decades ago in a slightly different 
system - the puls ed p arametric oscillator with a Kerr nonlinearity, where the classical and quantum dynamics were 
compared [^.IdOlli^. 

The field state becomes squeezed in the p^-quadrature for initial times, but the squeezing disappears for larger 
times [15, contrary to the ideal DCE case when (Apr)^ decreases exponentially with time [15• In the presence 
of damping (shown by the line 5 in figure [5]) the photon generation is still possible, but the oscillations of {hr), Q 
and (Apr)^, including the collapse-revival behavior, disappear [45- Moreover, the asymptotic value of the Q"fnctor 
(shown by the dash-dotted line in figure [5] and that can be calculated exactly [15) differs substantially from its value 













800 



FIG. 3: Time behavior of (fir), Q and (Apr)^ for the initial vacuum state, qr = 650|ar|, k = 0 and different values of uir- 
Notice the irregular collapse-revival behavior of (fir) and the maximization of the average number of created photons for 

(max) ry-i r\ 

OJr = —TlOar. 


during the transient, meaning that the field state for initial times is quite different from the asymptotic state. We 
also investigated how the dynamics is modified if the initial state is slightly different from the vacuum state. This can 
occur in actual experiments at finite temperature, so we considered the initial thermal state with the average photon 
number n, described by the density operator p = Pn\n){n\, pn = n^l{n + 1)"+^. The dashed line in figure^] 

shows the dynamics for n = 0.1 in the absence of damping, which should be compared with the line 3 calculated for 
the initial vacuum state. We see that for initial times the differences are very small and the oscillations of quantities 
(fir.), Q and (Apr-)^ persist. Therefore minor deviations of the initial state from the vacuum do not pose a serious 
threat on the experimental verification of the nonlinear DCE. 

From figures O and [3] we see that certain values of uJr maximize {fir) for times qrt < 5. We denote this value by 
^(max)^ noting that for another time interval the value of uir that maximizes the average photon number may be 
different. Since in actual implementations it might be difficult to maintain external modulations for a long period of 
time, the choice qrt < 5 seems appropriate to reflect the experimental reality. Figure |4^ shows the behavior of 
as function of qr/\ar\'. for large values of qr/\ar\ it is roughly given by = l.OQfjr- The maximum number of 

photons {nr)ra_ax '"^hen the effective detuning is adjusted to is shown in figure | 3 } d : for qr/\ar\ ^ 1 it grows 
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(d) 


<d~ 

V 

d 



FIG. 4: a) Behavior of the effective frequency (that maximizes (nr)) as function of qr/\ar\- b) Behavior of as 

function of qr/\Q.r \ for the effective frequencies and = 0. c) Behavior of (Apr)^;^ for these effective frequencies, 

d) Behavior of Q/ (nr) at different time instants. Curve 1 (2): tUr = (wr^'^^^) and the time instant of the minimum 

value (Apr)^. Curve 3 (4): ujr = (wr^'^^^) and the time instant of the maximum value of (hr). 


as (nr)jjjax = 1.36gr/|ctr|- On the other hand, if we ignore this optimization and set the value of Wr to the 

standard DCE resonance (without the Kerr nonlinearity), w) '^ = 0, then (nr)max grows linearly but with a 

smaller slope: (fir)max = 1.09(7r/|ctr|. Hence for large values of qr/\ar\ the optimization can significantly 

enhance the photon generation, facilitating the experimental verification. The downside of optimizing the effective 
detuning to increase (nr)max squeezing is decreased. FigureHJ: shows the smallest value of (Apr)^ achieved 

for Qrt < 5 when Wr is set to or Although in both cases the Pr-quadrature becomes squeezed, for 

^(DCE) squeezing is significantly stronger than for . 

Finally, in figure HJi we study the cavity field statistics at different time instants by plotting the quantity Q/{'hr) 
that quantifies the spread of the photon number distribution. We recall that Q/ifir) = 0 for the coherent state, 
Q/{nr) = 1 for the thermal state and Q/ifir) = 2 + l/(nr) for the squeezed vacuum state. The states for which 
Q/{fir) > 2 + 1/{fir), called “hyper-Poissonian” in [b^, have photon number distributions distinguished by very long 
tales with low probabilities that nonetheless cannot be neglected, so the mean photon number does not characterizes 
well the total distribution Q. The curves 1 and 2 denote the value of Q/{fir) at the time instant of minimum (Ap^)^ 
(shown in figure |3t) for and , respectively. We see that for we have approximately the squeezed 

vacuum state with Q/{nr) ~ 2 + l/{fir), while for we obtain a hyper-Poissonian state with a rather broad 

photon number distribution. Lines 3 and 4 denote the value of Q/ {hr) at the time instant of maximum {hr) (shown 
in figure HJd) for and respectively. In this case the field states are not very different one from another 

and have a super-Poissonian photon number distribution with 0 < Q/ {hr) < 1. 

Summarizing, in the presence of the Kerr nonlinearity one can optimize the photon generation from vacuum by 
adjusting the resonance shift (directly related to the effective detuning uir) as function of a,., and this is the second 
reason for the introduction of C, in equation (1131) . On one hand, this optimization decreases the amount of squeezing 
[^ . but on the other hand it can be used to produce novel field states @ very different from the squeezed vacuum 
state produced in standard DCE. 
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B. External classical pumping 


If the cavity undergoes a classical pumping one must add the term id (d^ — d) to the Hamiltonian ([8]). In terms of 
the auxiliary operators A and B we should add the term 


Hd = - 




1 

90 


-it(A-+2S-)/2 I 9o it(A-+2S-')/2 

+ ^2 e 


2A. 


^U)^-^t{A+/2-S+-r,(^'>) g-it(A.+2,5_)/2 _ git(A^+25_)/2 


B + h.c. 


(34) 


to the RHS of Hamiltonian (fT4l) . To the lowest order in go/A-, for the pump frequency = (wq + “ <^+ “ ~ C) 

one can simply add the effective pump Hamiltonian 


Hp — i 


1/ 


ip)* 




_^(p)g-. 0 ,/ 2 d^ 


(35) 


to the RHS of equation (|25]) . where the operator dr was defined in equation (l24ll . 

Neglecting the Kerr nonlinearity, the optimum resonance shift for DCE is C = 0. For the simultaneous modulation 
of the system parameters and the external pumping we obtain the general Hamiltonian of the form 

dApump = Qd,r —^dr “t“ h.C. , (1^6) 


where we introduced arbitrary time-independent complex coefficients g and In the Heisenberg picture the solution 
for the Hamiltonian (1361) is straightforward: 


dr = dr(0) cosh(|^| t) - i^d);(0) sinh(|^| t) -i- |[cosh(|^| t) - I] - sinh(|^| t) . 


(37) 


For the initial vacuum state, dr|0) = 0, we obtain for the average photon number 

{fir) = sinh^d^l t) + 2||L [coshd^l t) - sinh(|^| t) sin {2(j)g - (()^)] [cosh(|^| t) - 1] , (38) 

where we defined the phases as g = |^?| exp(j0g) and ^ = |^| exp(z0j) . 

For initial times, \^\t <C 1, we obtain 


{nr)^m^ + \9\^)t\ 


(39) 


so the average photon number does not depend on the phases (j)g and (f>^. However, for larger times the phases become 
very important as one gets 


(n) = sinh^dCI t) + 2||Le=^l«l‘ [cosh(|^| t) - 1], 

where k is any integer number. In particular, for large times |^|t ^ 

I , for 2(j}g 

+ for2</>,- 


for 2(j)g -41^= + 27 rfc , 

1 we get 

— (j)^ = —f -I- 27rfc 
4>i = ^ + 27rfc 


(40) 


(41) 


Therefore, by adjusting carefully the phase difference (2(/)g — (j)^) one can significantly amplify the photon generation 
with the help of one-photon pumping. Moreover, one could verify our model experimentally by measuring the 
dependence of {fir) on the phase either of the pump or the modulation parameters defined in equation ([5|). 


IV. DCE-LIKE BEHAVIORS WITH ATOMIC CLOUDS 


Besides forming the base of the toy model for Cavity DCE, for TV >> I the Hamiltonian (j4]) also describes accurately 
the interaction between a cold “atomic cloud” (e.g., polar molecules [12, US]) or ensemble of superconducting qubits 
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and a single cavity mode. In this case all the parameters of Hamiltonian ([8]) are controllable, and novel regimes of 
light-matter interaction can be implemented by modulating the system externally according to the law of motion ([5]) . 
The full solution is given in|Al and in this section we highlight the regimes in which excitations can be generated from 
vacuum. In section FVB I we shall describe another regime when pair of excitations can be coherently annihilated due 
to external modulation, in what we call “Anti-DCE”. 

In the dispersive regime, for the modulation frequency 

= 2(Ho - - C) , (42) 


we obtain the effective Hamiltonian 


Heff = 




5^ - 2 e 


iti. 


go 

A_ 


go 


+ ) +h.c. 


+ H]\[c , 


(43) 


where the non-Gaussian part ffjvG is given by equations (IA29I) and (IA34I) and we defined the time-independent 
parameter 


qi = - 


L 

T 


2 X 0 f 

A—y A..|. A.|- y 2Ho 



(44) 


Neglecting the non-linear terms and considering |i5x| <C |(?/|, we can write A{t) ~ A(0) — e A-)[B{t) — .B(O)], 

so to the lowest order in go/A- we obtain the total effective Hamiltonian (for the initial zero-excitation state) 


Heff ~ -f h.c.) -f d_[(HiH)2 - 5^5] . 


(45) 


Hamiltonian (HSl) is analogous to the DCE Hamiltonian but with the matter operator B instead of the cavity 
operator A. So this behavior corresponds to the DCE with matter, when pairs of atomic internal excitations are 
created from vacuum instead of photons (recall that the CM motion of atoms is prescribed externally). Notice that 
there is analogous Kerr nonlinearity term 6-{B^B)^, yet for \qi\ ^ |(5_| many matter excitations can be created 
from vacuum. We call this behavior ''''Inverse dynamical Casimir effecV (IDCE), since figuratively this phenomenon 
corresponds to exciting the internal degrees of freedom of the moving dielectric slab instead of creating photons. In 
section rV HI we shall also describe the “Anti-IDCE” phenomenon - an analog of Anti-DCE for the atomic degrees of 
freedom. 

Eor some modulation frequencies one can achieve simultaneous excitation of the cavity and the atoms. In the 
dispersive regime this occurs for the modulation frequency = A+ — 2(5+ — — C, when the effective Hamiltonian 

reads 


^eff = 


^itSy ffO ^;|2 _ ^—itSy 90 

' A_ A 


(ZMe-*‘^ I -AB + ^B'^ ) + h.c. 


+ ffjVG(go) + HNaiSg) ■ 


(46) 


The non-Gaussian parts are given by equations (IA29I) - (IA35I) and 


go 

qM = -^ 



AM) 

CUJ 

JM) 

1 

AM) 

AM)! 

[1* aJ 

A+ 

' “ A+ 

go 

+ A_ 


(47) 


We call this behavior “Mixed behavior”, since the photons and atomic excitations are created at the same rate 

\qM\ > |g/|- 

In the resonant regime, A_ = 0, excitations are generated from vacuum for the modulation frequency 77 ^^^ = 
2uJo — 2(5+ — (5^ — The total effective Hamiltonian is 

Heff = goe”*“^ _ QiR)^2itgo'^ 2.B 

+ 0We-2it9o 0(«)g2»tso^ ^2 

+ (- 00 ^^ + B^ + h.c. 


( 48 ) 
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The time-independent coefficients 0^^^ are given by equation (IA25I) and we neglected the non-Gaussian terms pro¬ 
portional to Sg. We see that for Q = 0,±2go one can create equal amounts of cavity and matter excitations. The 
photon generation for C = 0 unknown until a few years ago [s^ [s^ , yet it appears naturally in our formalism, 
as well as the non-Gaussian terms on the last line of equation (1481) . The detailed analysis of Hamiltonian (|S]) in the 
resonant regime (without the non-linear terms) was studied in |55l. in an attempt to describe the detection of 

DCE using small induction loops modeled as LG contours. 


V. NONSTATIONARY CIRCUIT QED WITH A SINGLE QUBIT 

Now we consider the limiting case A = 1 to study which phenomena exist for the most basic type of light-matter 
interaction under nonstationary conditions. From the practical point of view this analysis is relevant because it 
describes actual implementations in the circuit QED architecture, where the parameters of the cavity and the qubit can 
be modulated in situ by external biases and the one-photon classical pump is implemented in a straightforward manner 
[§-[l3,lll|. Nonstationary circuit QED has been studied in numerous papers duriiig the last decade (sl. 1^. . 

but here we generalize the previous results by working in the dressed-states basis |28l | and considering the simultaneous 
multi-tone modulation of all the system parameters. Moreover, we predict the new effect in which pair of excitations 
can be coherently annihilated due to external modulation, in what we call “Anti-DCE” behavior. 

As shown in |B] the wavefunction corresponding to the Hamiltonian (|3]) can be written approximately as 

CO 

IV’(Q) e“*‘^“6o(t)|(po) + X! X! ■ (49) 

Here |‘/5n,s) and Xn,s are the n-excitations eigenstates (also known as dressed states) and the “corrected” eigenvalues 
of the bare Jaynes-Cummings Hamiltonian 

Hjc = Won -I- Ho|e)(e| -I- go{da+ + a’l’d.) . (50) 

Coefficients b represent approximately the probability amplitudes of the dressed states and the index S labels the 
different eigenstates with the same number of excitations. The corrected eigenfrequencies and the eigenstates read 
approximately 


Ao ~ 0 , A„>o, 5 ^ Won- — +S-Pn , Pn = ^-+ ^gln , 5 = ± 


Iv’o) = \g, 0 ) , |</?n>o,s) = n) -b c„, 5 |e, n - 1) , 

where A± = wo ± Ho and we introduced the notation 

— cosOjYi^ Cyj.j^- 1 - — cosOjYi^ — sin( 


0m>o = arctan ■ 


A_ 


2goy/rn 


(51) 

(52) 


(53) 

(54) 


A. DCE behavior 


For a single modulation frequency matching the DCE resonance, ~ 2a;o, the probability amplitudes obey the 
differential equation [see equations (IB19I) - (IB25I) ] 


bm.T— y~![Q 








'm+2,S 


m.S.T^ ^rn-2,S\ ) 


(55) 


where we use the shorthand notation 6 o,t = ^o, -^o.r = -^o; 0 !i^ 2 T 5 ^ time-independent coefficient given by 

equations (IB20I) and (IB21|) . In the argument of the exponential functions there is an intrinsic uncertainty we call 
“Systematic-error frequency shift” (SEFS) due to the involved approximations. The estimative of is given 
in m since its order of magnitude is important to tune precisely the resonant modulation frequency. 
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The frequency matches the difference Xm+ 2 ,s — Xm,T only in the dispersive regime, |A_|/2 ^ go^/n, where n 
is an integer representing the number of excitations. Introducing the “detuning symbol” T) = A_/|A_| = ± we can 
write the corrected eigenvalues and eigenstates as (we denote Aq = Ao,i)) 

d s 

Am, I? = Wgm - am^ -5+ -^ , Am>o,-x> = Wgm - A_ + arn^ - 5+ + -A (56) 

Wm,^) ^ |g, m) + |e, W - 1) , Wm>Q-v) ^ -P(|e, w - 1) - |g, m)) , (57) 

where the effective cavity frequency ujg or We and the intrinsic “frequency shifts” are 

ujg=u:o + 5- - 5+ - 5^ , We = wo - <5- + <5+ - , <5± = ^ (58) 

A± A+ 

Hence the modulation frequency can couple either the dressed states |(/3m,xi) ^ \'Pra± 2 ,v) (where Iv^o.i?) = It’d)) 
or |i^m>o,-x>) \'-Pm± 2 -v)■ The former case occurs when the atom is predominantly in the ground state, so we call 

it g-DCE behavior. The second case corresponds to the atom predominantly in the excited state, so we call it e-DCE 
behavior. 

1) g-DCE behavior. Under the approximations and SETS 

\^m+ 2 ,-Ti,T>\^ \^m+ 2 ,v,-v\ |A-|, | 0 m+ 2 ,-X>,-X)l l*^-! (^9) 


l®m+2,-X),X>l^/|A-|, \^m+2,V,-v\‘^/\^-\> I ®m+2,-X>,-X) I ^ /1 A- I 

we define the effective probability amplitudes Cm as 

bo , m = 0 


\ bm,v , m > 0 

Adjusting the modulation frequency to 

77^^) = 2(a;g - C) , 

where C is an adjustable resonance shift, we obtain the differential equations 

Cm = |d+| + l)(m + 2)e-2*‘[f-2“(-+i)lcm+2 

-e-*^+ v'TO(m - l)e"**[^-2“('"-i)lcm-2 


(60) 


(61) 


(62) 


(63) 


= - 



■ Xo^ 

\ ^ 1 A 


2A Ho 

.(D) 1 

^9 1 ^JD) 

u+ 


\ 0 — 10- 
' Wo 

A+ 

A+ 

-1 

0 

h 

o 


= |d+|e*^+ 


(64) 


Comparing with equations (|2^ and (|2^ we see that under a trivial phase rotation dr —)• the dynamics 

of Cm is described by the nonlinear DCE Hamiltonian (1251) with ujr = (, ar = —a and qr = |i?+|. If we also apply a 
one-photon pump with the frequency = ujg — then one can simply add the term — h.c.] 

to the RHS of equation (1251) under the additional approximations 


£d'/^ £dgo 


0(<'r) 



2) e-DCE behavior. On the other hand, for the modulation frequency 


77^^) = 2(we - C) 


(65) 


(66) 
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we define Cm = bm+i.-v and obtain the differential equations 


= |0_| v'(m + l)(m + 2)e-2*‘[^+"“('"+2)lc™+2 


(67) 


where 





2(5_»o 

A+ 50 


+ ie 


(D) 




( 68 ) 


These results are valid under the approximations (l5^ - (l60)) with replacement V —>■ —V. So the dynamics of Cm is 
again described by the nonlinear DCE Hamiltonian with Wr = C + 2a, a^. = a and qr = |'(?-|, and for the external 
pump with frequency = uje — C one can simply add the term —+ h.c.]. 

Thus the nonlinear dynamical Casimir effect exists even for a singe qubit, so it is an intrinsic phenomenon of the 
light-matter interaction in nonstationary systems and can be observed in the circuit QED architecture. For a single 
qubit there are two possible modulation frequencies, equations (1621) and (1661) . whereas for N ^ 1 we found only one 
resonant modulation frequency, equation (US. The origin of this apparent discrepancy is trivial: in section IIIII we 
assumed that {b^b) <C N, so the case when the atoms were initially in the excited states was automatically excluded 
from the treatment. The photon generation rates |i?+| and |i?_| are of the same order of magnitude, but may differ 
due to the phases for simultaneous modulation of several parameter. The effective detuning tOr and the Kerr 
coefficient ar depend on the initial atomic states, so the resonance shift must be adjusted accordingly to optimize the 
photon generation. By increasing the number of qubits we simply make the replacements |a| —>■ fV|a| and 6- —>■ NS-, 
as can be seen from equations (fT51) . (IMl) and (1551) . So for = 0 the maximum number of photons created from 
the initial zero-excitation state (shown in figure |3 }d) is not altered by increasing the number of atoms, although the 
photon generation rate undergoes a TV-fold increase. 


B. Anti-DCE behavior 


In the dispersive regime, for the modulation frequency (we neglect the Kerr nonlinearity a to simplify the expres¬ 
sions) 


5m^ — 2wo T — 3(5^ -|- 2 (5_ — (5_|_) {M 1) , 


(69) 


where M is a positive integer, we obtain the differential equations for m > 0 


i _ 


— (5_|_) (m—M) ^ 


'm+2,'P 


^m+2,'D 


~ ^m+2,-V,V^ Um..- 




(70) 


0^j2,-x>,x> = ^ 2a;lA^° + ^)(^ + ^) 




-k 


Wo + A ^ 




2^0 “ I ” 2^0 “ 1 “ ^0 .^0 


(71) 


The involved approximations are 


l^m-|-2,-D,-r>h l'^m+2,'D.'Dl> \^m+2,T>,-V 


«|A-| 


(72) 


0(^i",V) - • (73) 

Under realistic conditions we have |(5_| ^ |0m+2 -v -dI’ amplitudes bM,-v and &m-i- 2 ,x> are effectively 

coupled. Therefore this modulation roughly couples the states |e, M — 1) \g, M + 2). In other words, for the initial 

state \g) ® couple the subsets \g,M) o |e, M — 3), thereby annihilating three photons (two 

excitations in total) via external modulation. However the coupling rate |0|(f+2 -v v\ very small, so the frequency 
(1691) must be fine tuned (taking into account the Kerr nonlinearity and SEES) and the transfer of populations between 
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the states takes a long time. If only ^ 0, as in parametric down-conversion, then 0^+2 -x> d — ^ process 

does not occur at all. Noticing that in the dispersive regime the coupling \g, m) O |e, m — 2) via one-photon pumping 
is prohibited, we conclude that the subtraction of system excitations via external modulation (when the atom starts 
in the ground state) only occurs for the time-modulation of parameters w, 17 or g. 

This phenomenon persists in the macroscopic case for a cold atomic cloud. As shown in lA 21 for the modulation 
frequency ~ 2ujq + A- = Swq —flo we obtain effective Hamiltonian of the form Heff — /{ANA-)]{A^B'^ — 

A^^B) + • • •, where other nonlinear terms are given by equations (IA29I) and (IA36I) and we neglected the contributions 
of and en. So when the atoms start in the ground states there is an annihilation of three photons accompanied 
by generation of one collective atomic excitation. Since in this case the photons are annihilated by virtue of external 
modulation of the system parameters, including the prescribed motion of the atomic cloud, we call this effect “Anti- 
DCE”. This name should not be taken too literally because such behavior cannot be implemented with a dielectric 
slab for which the parameters H and g, and hence the resonant modulation frequency r]^^\ are not known. 

In section Irvl we described the IDCE behavior, when pairs of atomic excitations are generated from vacuum for the 
modulation frequency ~ 2ilo. By symmetry in Hamiltonian ([ 8 |), there is also the “Anti-IDCE” behavior, when 
three atomic excitations are annihilated (for the cavity field in the ground state) due to the modulation of system 
parameters with frequency pe 2Ho — A_ = 3Ho — wo- This effect is described by the effective Hamiltonian 

Heff — —[i5-£^^^'’/{ANA-)]{I3^A^ — B^^A) -\ -, as given by equations (jA29l) and (IA37I) . In practice the Anti-DCE 

and Anti-IDCE behaviors are very difficult to observe because the involved coupling rates are quite small. However, 
they are interesting from the purely theoretical point of view for constituting examples of motion-induced coherent 
annihilation of excitations. 


C. Generation of entangled states 


Now we briefly review some practical schemes to generate entangled states in circuit QED with time-modulated 
parameters, studied previously in [H, 113,165. In the dispersive regime, for the modulation frequency (neglecting the 
nonlinearity a) 


=A+-2{S--S+){M + 1)-S^ 


and approximations 


|0 


(S) 

m+ 2 ,'D,’D 


’ l®m+2,--D,-E>l> l®m+2,-X',-X'l l^+l 


- 1 © 


m,T 


m+2,X>,X> I 


VlA+i, |0!;f|2,_p,^l7lA+|, |0L^|2 .-i,.-i,IV|a+| 


m+2,—X>,— 


we obtain the equations (denoting 60,15 = ^ 0 ) 


:>(S) 


6m,15 — 0m+2,I5,-I5® 


l>m+2,-V 


-75 = —0 


^m+ 2 ,—X> 

it2{S--S+)(M-m)i 
m+2,V-V^ <5m,75 


aiS)* 


(74) 

(75) 

(76) 


(77) 


0 


(S) 

m+ 2 ,I>,-'D 


-75 = -^9oH\/m+ 1 


/i7i/ \ 

\ ^ _£3_ 


(S) 


SS) 


- 1 + 


2 e 


iS) 


A_ / A 


90 


(78) 


/ C\ 

For |6_| ^ |0m7 75 - 75 I only the amplitudes 6 m ,75 and 6 m+ 2,-75 are effectively coupled, so this frequency roughly 
couples the states |g, M) o \e,M + I). For this reason such behavior was called “AJC regime” in [13, [EHl and “blue- 
sideband transition” in (53| . recalling that the Anti Jaynes-Cummings (AJC) Hamiltonian is Hajc oc ad- -|- a^(t+. 
For N fA A this behavior turns into the “mixed behavior” described approximately by the effective Hamiltonian 
Heff oc AB + AAB\ as follows from equation (1461) . Moreover, for the external one-photon pumping with frequency 
= Do — ( 6 - — 6 +) (2M -I- 1) (neglecting the nonlinearity a) we obtain 

(P) 

h - -2it{S_-S+){M-m)i 

^m,T> — ^ 2 |A 

(p)* 

i) .T r, - i^Ald _ 

















16 


TABLE I: Abbreviation of the effects with atomic clouds in the dispersive regime. The asterisk marks the effects that persist 
for a single qubit, stands for the approximate modulation frequency. 


Abbreviation 

^0) 

Main effect 

DCE (*) 

2 ujq 

Generation of pairs of photons 

IDCE 

2Do 

Generation of pairs of atomic excitations 

Mixed (*) 

Wo + ^0 

Equal generation of photons and atomic excitations 

Anti-DCE (*) 

3wo — 

Annihilation of three photons 

Anti-IDCE 

sr^o — ^0 

Annihilation of three atomic excitations 


under the approximations ed^/m |A_| and 0{v^j-) ^ e^/|A^|. So for |(5_| 5 oed/|A_| one couples only the 

amplitudes 6 m, u bM+i,-v, corresponding to the selective excitation of the atom \g,M) o |e, M) conditioned on 
the presence of M photons in the cavity field. From the first line of equation (IB19|) we see that in the dispersive 
regime one can also couple the amplitudes bm,T> and bm,-v, or roughly the states with the same number of excitations 
\g,M) o |e,M — 1), by employing the modulation frequency ~ A_. This behavior was called “JC regime” in [^1^ 
and “red-sideband transition” in |54| . The generation of a single photon from vacuum and the transfer of populations 
between the cavity field and the atom using the red- and blue-sideband transitions was studied in details in . 

In the resonant regime, A_ = 0, we can couple the dressed states ^ \g^m+ 2 ,s), where \g}m,±) = 

{\g,m) ± |e,TO — 1)) /\/2, for any values of to, T and S by the modulation frequencies = \m+ 2 ,s — The 

corresponding coupling rates are of the same order of magnitude for any T and S and are given in IB II Besides, the 
states ^ \^m+i,s) can be coupled by the classical pumping with frequency = Am-i-i,5 ~ ^m,T- Therefore 

combining the temporal modulation of the system parameters with the external one-photon pumping one can create 
arbitrary superpositions of dressed states with a high degree of control. Moreover, one could apply several resonant 
modulation frequencies at once to study the dynamics under the multi-tone modulation, when many dressed states 
are coupled simultaneously with controllable rates. 


VI. CONCLUSIONS 

We showed analytically that Cavity dynamical Casimir effect is contained implicitly in the most basic form of 
the light-matter interaction - the dipole interaction between a single atom and a cavity field mode under external 
modulation of the atomic parameters. This phenomenon is intrinsically nonlinear due to the nonharmonic energy 
spectrum of the atom-field system, so the number of photons created from vacuum is limited and the resulting field 
state can be quite different from the squeezed vacuum state. The atom becomes entangled with the field and the 
average photon number exhibits collapse-revival behavior as function of time, very sensitive to small shifts in the 
modulation frequency. The effect persists when the number of noninteracting atoms N is increased, and for N ^ 1 
our approach behaves as a toy model for a oscillating dielectric slab inside a stationary cavity. The precise knowledge 
of the atomic parameters is not required to achieve DCE, since the photon generation occurs for the modulation 
frequency in the vicinity of 2ujo. Moreover, for additional external classical pumping the photon production via DCE 
can be substantially enhanced for appropriately chosen phase of the pump. 

If the atomic parameters are known and controllable in situ our model describes the nonstationary circuit QED 
architecture (for N 1) or cold atomic clouds (for N ^ 1). In this case we can employ other modulation frequencies 
to realize new effective regimes of light-matter interaction. In the dispersive regime these frequencies and associated 
effects are summarized in table HI In a new effect, that we called “Anti-DCE”, the modulation of atomic parameters 
can lead to coherent annihilation of three photons accompanied by the creation of one atomic excitation; however, 
the associated transition rate is very small so this behavior hardly can be implemented experimentally. Besides, we 
found atomic analogs of the DCE and Anti-DCE behaviors, when the photonic and the collective atomic operators 
are interchanged. Einally, we demonstrated that entangled states (dressed atom-field states) can be generated in a 
straightforward manner in nonstationary circuit QED. 
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Appendix A: Analytical results for A ^ 1 


Following the method described in [2^ we write the solution for the annihilation operators a and b in the Heisenberg 
picture as 




a = 


{(i+cLh + gobh)e + {P-clh - gobh)e'*-^/'^ 


(Al) 


b = 


„-ztA+/2 




iP-bh + goah)e + {l3+bh - 


where we defined time-independent parameters 


/3 = ^A2_+4g2 , A± =cco±Ho , /3± = ± 


(A2) 


(A3) 


so that A_ stands for the bare atom-cavity detuning. 

The auxiliary annihilation operators at and bh satisfy the bosonic commutation relations [a;i,aj)] = 1, [bh,bl] = 1, 
[ah,bfi] = = 0. Under the approximations 


£uj £n A_ Ed (S’l'S) 


11 < I 

go ljq luo luo luo ^ P 


they read 


LUq LUq OJq 




(A4) 




(A5) 


Here E = (a^a) + {b^b) denotes the total number of excitations in the atoms-field system and we defined small 
“intrinsic frequency shifts” 

A = ~ 

A+ ’ ' “ /?" 


h = ^ = ^ , Ss = Y.^ - A_4^)) . (A6) 


6- is the standard dispersive shift, is the collective Bloch-Siegert shift, S-^ is the shift due to the term ixo{^^^ — 
and 6s is the shift due to a possible modulation of the system parameters with a low modulation frequency « p. 

The independent annihilation operators A and B, that also satisfy the bosonic commutation relations, are defined 
implicitly in terms of the small time-dependent functions I-U4 I, |-Fb|, |J^_ab| 1 


= E 


/ 1 

j w 


(£^0^ - 2e^l'>~gl + ^ffoA_) 




g 


( 3 ) 


pit{P^’+P) — \ ... ... -P) — I 

+ P- 9 I - Sg^gO^-) - , o - + (^- 50 - 4^^5oA_)- 73 -^- -b C.C. 


gU) + p 


gU) — p 


(A7) 


= Ej^ (en4^+2e45o - 24^^A_5 o) 




g 


( 3 ) 


-p-gl - 4^4oA_)- . . -(£^4o - 4^ 5oA_)- 7 -^-^— + C.C. 


g^A + p 


gU) — p 


(AS) 


^ab = E 


' 1 


(e(f4oA_ + 44^'4 o) 




3 2/32 

+(e^4o - 4'’^A_) ( P- 


+ (e4*5oA_+4£4>52) 




gU) ^ g^A 




g(A — p 


-P+- 


gU) + p 


S 3 )*~ -AO*. - 1 ^ - 1' 


+ (e_ *50 - 4'’^*A-) P-- 


g(A -\- p 


-p+- 


gU) — p 


(A9) 
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The time evolution of the operators A and B is governed by the Heisenberg equation of motion idO/dt = [O, Heff], 
where O = A, B and the effective Hamiltonian can be written as 


^eff = Ha + HNcigo) + ^ Hno^sx) ■ 

X—g,uj,Cl 


(AlO) 


Hq denotes the Gaussian part containing linear and quadratic combinations of A and B. HNcido) and Hng{sx) 
denote the non-Gaussian parts, of the fourth order in operators A and B, proportional to go and ex, respectively. 
For simplicity we shall consider only the term Rxcisg) in the last term of equation (lAlOl) . since the resulting general 
expressions are too long to write out explicitly. 

Eliminating the rapidly oscillating terms via the Rotating Wave approximation (RWA) [2^ we obtain for the 
Gaussian part 




+/ 3 _ (e^l^*go - Ab^ 


--e 

2 


J2it5s 


-~goQo^ + + /3+0^'^e" 


i/3 


+e' 


it5-y 


B^ 

-iA_0[,^^ + 500+^-**^ - 5o0-^*‘^ inj 
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-it 


g ..[X+/2-S+-A^'') |gWgit(<5x-^o) 


/3_e 


itp/2 


- p+e 


-it/3/2 


~goe-^^f</^-~goA^P/^ 


ii| + h.c. 


The time-independent coefficients are 


0« = - 2xoV^^^ , 02'^ = 


A" 


(All) 


(A12) 


^0^'^ = ^ 


Si) 


Si)' 






90 


77(1) — (3“^ 


{es’go - 


(A13) 


yii) _ 

Vq — 


' Si) n 
gx go 

Xo 


' h pA) 

9o^+ 

nii) 


V 


ii) _ A_ 


jjU) — p2 


{e_go - 


= 


g-’» - T (4R. ± 24>>9» ± 4> A-)^ 


±£q^90 


(A14) 


(A15) 


vA'’ = — 

± 732 


^‘'xP± 

Xo 


^ ^(eL^'^/3± + ± 2e2^ffo)^ - (e^'^A. - 


(A16) 


Due to elimination of the rapidly oscillating terms one introduces intrinsic uncertainty Av in the arguments of the 
exponential functions in equation (lAllI) of the order 


0(At/) 


I1^9L(£\ XA(l. 

p) iOo\P) Wo’ Wo’ Wo V'^o/ ’wo’woVwo 


(A17) 


We call these contributions “Systematic-error frequency shifts” (SEFS), since they appear due to systematic simplifi¬ 
cation of the differential equations for dh and 6/j. In practice SEFS slightly alter the resonant modulation frequencies 
that give rise to nontrivial behavior, so ultimately they must be found experimentally or numerically. 
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Under the additional approximations 


eVe 

N 


<C 1 (Resonant regime: A_ = 0) 


go Ve 


^ ^(ata)(6tS), ^(ata)| < 1 (Dispersive regime: |A_| /2 » go) 

one obtains for the first non-Gaussian term in equation (lAlOl) 

HNGigo) = -^|| [3g2A_it2i2 ^ 2A_ (A^ - 2g^o) AB^B 

+8go (A?_ + go) - G^qA- 


(A18) 


(A19) 


(A20) 


-Ago {A^_ - 2g^) it2^5g-ii(5x-2U) _ ^2^ ^12^2 


+ h.c. 


The simplified expression for the non-Gaussian term Hj^Gi^g) strongly depends on the modulation frequency ^ as 
can be seen from equation (|S]), so we do not write it explicitly due to its length. In IA 21 we shall give the approximate 
results for H^fcisg) in the dispersive regime for high modulation frequencies g^A ^ 2a;o. 


1. Simplified formulae in the resonant regime 

For A_ = 0 we obtain the simplified expressions 

a ~ cosgot — singot 

I, ^ ^-^t(uo-s+) [Be**^''=cosgot-*ie*‘(‘^^-^'’)sin5ot 


The Hamiltonians are 


Hn ~ 


g-it{2go+2Ss-Sx-v‘-^'‘) _|_ g(l)* g** (2go-2U-l-<5x^ 
_1 g-it(wo-i5+-J)'^’) |gWgit(<5x-'5s) ^g»*So _|_ g-**So^ ^ 

+ ^g-iiso _ ^tgo) 

_j_^^ g-it(2i^o-25+-r;<^>) Igit^x — 0^A ^it2go'j 

_|_lg2®‘(‘5x-'5s) -I- ©t()g-**23o _|_ 0(l)git23o^ ^2 

+ ig2iU. |^_©(f) + ©Cj)g-*t2go 0U)git2goJ ^2 | _^ ^ 


1 AB 


HNcigo) = [a'iA + B^B - l) ^t5g-*t(5x-25.) 


The time-independent coefficients become 


©(U = P— (_ pU) 
2~goWA^^ ^ 


■ 9o 0) 

— » . . -n'o £- 


©« = i 


(?7(U)2^- 

1 f^gO-)±,gO-A £ 
giA \go + J 


(A + „■ g(l) 
X ^ ^^9 


90 


(A21) 

(A22) 


(A23) 


(A24) 


(A25) 
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2. Simplified formulae in the dispersive regime 


In the dispersive regime, |A_|/2 go, we have (3 ~ |A_| + 2|(5_|, where S- = is the collective dispersive 

shift. The operators read approximately 


g ^ ^-it{ujo+Ss-S+-S,i^) 


A - sin (A_t/2) 


Introducing the “detuning symbol” 


B - sin (A_t/2) 


=± 


(A26) 


(A27) 


the Hamiltonians can be written as 


Ha = \vY, 


A_ 


0) go _ Hi) 1 -it(A_+2A+25.-5,,-Pr,0)) 


e 

V 


.h. (_ £(f) 
A2 A_ 9 


^it( A_ +25_ -25„+5,, ) 


-V 


ini 


,1 V' 

9. 


^U)^it{S,^-Sa) I g-it(A_+25_)/2 _^gZt(A_+25_)/2 j ^ 


+ ^g-it(A_+25_)/2 _ git(A_+2A)/2^ ^ g-zt(A+-25+-,,0)) 


+Pe^ 

1 


itS-^ 


2it5s 


-i|A_|0“ +5o0^'^e-**(^-+2^-) -go0^],e**(^-+2^-) 


A2 
AB 

\ A h.c. 


(A28) 


HNcigo) — 


2N 


B^B^ _ _ 3^^t2^2 _ 8^it5t^2g-z4(5^-25.) 

A^, A_ 


I 0 9o j^2^-2it{5,,-2SA _|_ 4_^^t2^^g-*t('5x-25„) 

A^ A_ 


+ h.c. 


(A29) 


In the first line of equation (IA28I) we introduced the notation (0)_,_ = O and (0)_ = O*. The time-independent 
coefficients in Hq are 


(a(A ^ 50 

° ~ A_ 


, 0 ) Jj) = 0 ') Jj)' 

j__ I AA + j££ 2_ _ j£ff_ + 2^ 

A_ / niA ■nU) Qo A_ 


(A30) 
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U) 

V 
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2xo + i5. 
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U) _ 2A_ 


- \ e. 


0) - A. 
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U) 

LJ 

lU) 


M) 


iS. 


u) - A, 


-iS. 


nU) _ 2A_ 


V 


V 


O') - A_ 


a _ Aj) 


90 


(A31) 


fiO) ^ -Ts 5o 




,0) 


-b f 


r^O) + 2 A_ £q i . yyO) -t- 2A_ 


— i- 


A_ J 77O) -b A_ 77O) -b A_ ryO) 77O) -b A^ 50 


xO) 

-g I 

A_ 


(A32) 


Neglecting the rapidly oscillating terms under the approximations (IAI9I) we can obtain particular expressions for 
Hng (Cg) for concrete modulation frequencies: 
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for the DCE modulation frequency ~ 2a;o 

AD) 


HnG i^g) = i 


c 

■ ^9 _ -it[2(^uJo+S--S+-S^)-g(^^] 


4iV go 


+A4^AB'<B^e- 

A_ 


_4_S!o _ 2_^i3^tgit(<5x-45.) 

A_ A?_ A_ 

~2 *"3 ~2 

_j_^_5o_^t^3g—2it(5s _ 2) | 2 ^t^3g~2^^f<5Y~flsl 

A2 /\3 ^2 


+ h.c. 


for the IDCE modulation frequency 77*-'^^ ~ 2120 

Ai) 


HnG (Sg) = i 


AN go 


-B^B^ + A^Ab^ - 


~2 " ~3 

_}__5o_^t^3g2z7((5;^—2(5s) I 2 A^ ^3^——2^a) _j_ 2_^0_^3^tg3z7((5;^ —2(5s) 

A?_ A„ A?_ 


I 1 ^4gt^2^zt((5v-2<5s) _ 5 ^0 42 j^tgg2it(ay-2^3) 

A_ A2 


+ h.c. 


for the mixed modulation frequency g^^'> « A_| 


dM) 


A (~ \ • 90 ^g —zi[A j. —2(5+—i5v — 

[ + J 


iSt52 _ ^^ti52g-zt(5^-25,) 

A_ 


~2 ~3 ~ 

_^^t^3g-2it(,5^-2,5.) _ 3^^t^3git(,5^-25,) , 90 ^■\^3^-it(5^-2Ss) 

A?. A?. A_ 

- 24 ^A‘^B^BAA^^-^^A + -^^3^tg2it(5x-25.) 

A^ A_ A^ 


h.c. 


for the Anti-DCE modulation frequency 77^"^) Ri 2(+;o + A, 


r ~ A^) 

HNG{eg) = 7^^£|_e-i‘[2‘^°+^-+4'5--2-5+-35x-’7''"d 


3 Dt 4it(5s 




4A A_ go 

-JL ^2g-it(5^+25.) 3|o_ _ ^1^2^^ g 


— 2itS-y 


+ io_ (#5- SiU) 52g-zt(35^-25.) _ J^^t53g-4zt(5^-5.) 


+ h.c. 


for the Anti-IDCE modulation frequency ra 2r2o — A 

AAI) 


HnG {Sg) = * 


.S- go eg _^ it[2no-A_-45_-2(5++5x-r)<'^'d 


-i1'i33g4it-5(» 


4A A_ go 

+ (siU - B+b) + 3 (^Ab^B^ - 

+ -^ (^A^A - 3B^Bj 42gZi(35^-2,53) _|_ .^^3^tg4*t(<5x-50 


^22t5v 


h.c. 


(A33) 


(A34) 


(A35) 


(A36) 


(A37) 


Appendix B: Analytical results for N — 1 


For A = 1 we work in the Schrodinger picture and expand the wavefunction corresponding to the Hamiltonian (|4]) 
as [H 


\lj;{t)) = e Ao{t)\ipo) + X! X! ® , 

n=15=± 


(Bl) 
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where An.s and are the n-excitations eigenvalues and eigenstates of the bare Jaynes-Cummings Hamiltonian 

Hjc = Won -I- f^o|e)(e| -I- go{aa+ + a^a-) . (B2) 

The index S labels the different eigenstates with the same number of excitations n. 

The well known eigenfrequencies are 

Ao = 0 , A„>o, 5 = Won -k ^ , Pn = 4ffo^ , 5 = ± , (B3) 

where A_ = ojq — Hq is the bare detuning. The Jaynes-Cummings eigenstates, also known as dressed states, are 

I'Po) = \g, 0) , |<Pn>o,s) = Sn,s\g, n) + Cn,s\e, n - 1) , (B4) 

where we introduced the notation 


Sn,+ — sindyi, Sn,— — OOsOn, — OOS^yi, Cn,— — siu 


with 


dn>o = arctan ■ 


2goy/n 

We introduce new time-dependent probability amplitudes b{t) via the relations 

Ho(t) = 


(B5) 


(B6) 


(B7) 


A, 


.r(t) = ^ 


—UJ,Q,g 


—Am.-r-t-’)'-’') _ 

(Am,r — ^m,-T + {^m,T “ ^m,-T “ 


-e-rk - 




X exp 


^e: e 


n: 


k,j 

m,T,T 


k—uj,Q,g 


V 


U) 


cos 'j _ gQg ^U) 


(B8) 


where the sum runs over “high” modulation frequencies g^A) > and we defined the time-independent coefficients 

Bmfr.s = (B9) 


Bmfr.s = £nw^n\^m,T\e){e\ipm,s) 
BmV.5 = ^g'>J^gH‘Pm,T\{aK + 0,^a-)\ipm,s) ■ 


(BIO) 

(Bll) 


These quantities are calculated in a straightforward manner using the dressed states. 

In equations (IB7p - (IBSp we introduced small “intrinsic frequency shifts” [2^ due to the elimination of the rapidly 
rotating terms throughout the derivation: 


,(i) 
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—iE kr E 


|si,i 


4 ^ I Ai,s -I- g^A 


(B12) 
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1 Jj) 
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5=± 

^m,T + V^A Am+i_5 — Xm,T + V^A 


(B13) 
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where we defined 


,(i) 
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m>l,T ^ / y 


L(5) 

V 

r2 r2 

5^1,m,5,r ^l,m+l,r,S 

fd 

5=± 

^m,T — Am-1,5 + Am+1,5 — Am,r + 


(B14) 


(2) _ ‘-2,s5o + 2®2,5^0 (2) _ ^m+2,r,s5o + ■^2,m+2,T,S^0 


5 =± 




2,5 


5=± 


'^m+2,5 


(B15) 
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^2,r “ 


^2,5,rffo + ^2,2,5,rAo ^m+2,r,55o + ^2,m+2,r,5Xo 


^m,r 


A, 


m+2,5 ^m,T 


(B16) 


^m>2 


,r=E 


5=± 


^'m,S,T9o + ^E.5,rAo ^m+2,r,5go + -^E+2,T,5Ao 


^m,T Am—2,5 


Am+2,5 Am,T 


(B17) 


Am+2,r,5 = (‘/^m,rl®'5’-|7^m+2,5) > Lk,m+k,T,S = (V’m.T 1^^ l7^m+fc,5) • 


(B18) 


The new probability amplitudes obey the differential equations (to simplify the notation we denote 6o,T = ^o, 
It’o.t) = Iv^o) and Ao,r = -^o) 


5m,r = 




k—uj,Q,g 




U) 

m-\-2,'T,S 


^-zt(Am + 2,S-Am,r-’) 


0)'| 


^m+2,tS 


^m,S,T^ '°m-2,S 

+ 2 / , / , ® ' '^tjl,m+l,T,50m+l,5 

5 i 


e V ’ ^ 7-0m-l,5 

with time-independent coefficients (where m > 0) 
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(i) 

2,r,5 
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2 A.,^-A.,5 + r;0) ' - 2 l4^'A2,r,5 - *5W^2,2,r,5) 


7^=± Z=ci;,r2,5 


(B19) 


(B20) 


0(5) - - V V 

'^m+2,r,5 “ 9 Z^ 


7^=± Z=a;,n,p 


g0^m+2,T,n — iXoL2,m+2,T,n^l,j 


\ \ _i_ />o(7) ^m+2,7^,»S 

Am-|-2,'R — Am-l-2,5 + ' 
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X4>\ 


U) 


g0^m+2,TZ,S — ‘1X0^2,m+2,n,S T-,IJ 
Xm,T — Am,TC + 


(B21) 


2 ^”1-1-2,T,5 — ^L2,m+2,T,s'j ■ 


Equation (IB19I) was deduced under the following approximations [recall that (/) stands for “high” modulation fre¬ 
quencies 77(5 ') > cuq and k = uj,n, g] 


Yikji _ Yfkji 

55m,5,5 51^ _5 _5 
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771 , 0,-0 


n 


kji 
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k,jl 


mit2,5,5 ^m,—S,—S 
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kjl 
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k,jl 
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< 1 


(B22) 


{g0)£g}5^m-|-2,r,5 

Wo 


{X0;5x}^2,m,r,5 

Wo 


5dTx m-t-l,7",5 
Wo 


< 1. 


(B23) 
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U 


Notice that in equation (IB19I) the resonant modulation frequencies correspond to the difference between two 
corrected” eigenfrequencies defined as 


V,r 


= A 


■m,r + 


^(2) 


(B24) 


(we denote Ao.r = Ao). So the Jaynes-Cummings eigenfrequencies are corrected by the frequency shifts and 

V- equation (IB24I) we neglected the additional frequency shift due to the modulation depths Eg, e^j, £n 
and of the order 


o(<V) 


mez 


UJQ 


(n: 


kjf 

m±2,tS,»S 


- n: 


k,ji ' 

m,—S,—S t 


Wo 


Wo 


for k = g 


(B25) 


We call these neglected frequency shifts “Systematic-error frequency shifts” (SEES), since they appear due to the 
systematic simplification of the differential equations for bm,T using the RWA [^. The knowledge of SETS is important 
because they slightly alter the resonant modulation frequencies, so ultimately they must be found numerically or 
experimentally. 


1. Simplified formulae in the resonant regime 

For A_ = 0 we obtain the expressions 

1 




Si) 12 




(B26) 


_ ^ ^ I ^ 

Xm>o,s = uJom + Sgo^/m-{6++mSx) - 


U) 12 


4 2^ Wo + 


(B27) 


A± ^ A+ 


(B28) 


(|5, m) + 5|e, m - 1)) 


(B29) 


The coefficients are (for m > 0) 


c,0) _ 1/2 


-O') Jj) 


(go - 2txoV2) ^+g,^-(l + eO) + 


jyO) J 
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(B31) 
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\/to -b 1 
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•v/m + 1 


qU) _ _ 


U) gO') 

-[50 + 2ixo{VrnT2- - 9o^ + £g^ + i£^^^(yrrr +2 - Vto) 


(B34) 
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\/to + 1 


0) pO) 

[50 - 2ixo(\/TO + 2 - + 90^ - 4'^^ + ie^^H'ym + 2 - ^/m) 


(B35) 


Li.i,T,S = , ii,m+i,r,r = ^ (vWTT + -v/r^ , Li,m+i,r,-r = ^ (\/m + 1 - y/m) . (B36) 


From equations (|B22p . (|B23p and (|B25p we derive explicitly the underlying approximations and SEFS in the resonant 
regime 


£oj,£Q,go'/M, egVM, Sd^/M, XqM, e^M < wq 


00/3) ^ 4 4 ^4 

Wo ’ Wq ’ Wo ’ Wo 


(B37) 

(B38) 


2. Simplified formulae in the dispersive regime 

For |A_|/2 goyjn we obtain after expanding /3„ in equation (jB3l) 

Xn.v — — ctn^ , Xn-v — Won — A_ — 5-n + an^ 


(B39) 


|‘/5m.p) ^|g,TO) + ^v/n|e,m - 1)^ , _p) ~-P ^|e, m - 1) - ^ Vm|g, m)^ , 


(B40) 


where V is the “detuning symbol”, equation (IA27I) . and the effective Kerr nonlinearity strength is a = gt/lxt .The 
coefficients are (for m > 0) 
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(B44) 


0 


0 ) 

m+2,X>,X> 


5-yJ{m + 1 )(to + 2 ) 

_ 2 A- 

2 *Xo/ e4 

ryO) _ 2A_ e4 ^ 

2 

1 

<1 

1 

_ 1 

5- j ?y0) ryO) — A_ 

ryO) — A- go ' 


(B45) 
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— Vm + 1 , — -rr—r ; -^ 1 , 771 + 1 ,- 15 ,- 1 ? — , -^^ 1 , 777 + 1 ,- 15 , 1 ? ~ ^[( 50 /^-)^] ■ (B47) 
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The frequency shifts are 
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jy^ 


( 2 ) 

m>0,P 


= -(to + 1)<5+ - (^m+^'^6^ , i^f^>o,-i? = (m - 1)5+ “ ^x- 


(B49) 


From equations. (|B22ll . (IB23|1 and (IB25p we derive explicitly the underlying approximations and SEFS in the 
dispersive regime 


eu,jeQ.,go\/m,egym,edy/m, xom,e^m <C wq 
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